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Abstract 

In this paper we present a rather general phenomenological theory of tick-by-tick 
dynamics in financial markets. Many well-known aspects, such as the Levy scaling 
form, follow as particular cases of the theory. The theory fully takes into account the 
non-Mar kovian and non-local character of financial time series. Predictions on the 
long-time behaviour of the waiting-time probability density are presented. Finally, 
a general scaling form is given, based on the solution of the fractional diffusion 
equation. 
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1 Introduction 



The importance of random walks in finance has been known since the seminal 
work of Bachelier [1] which was completed at the end of the XlXth century, 
nearly a hundred years ago. The ideas of Bachelier were further carried out 
by Mandelbrot [2], who introduced the concept of Levy flights and stable 
distributions [3] in finance, and by the MIT school of Samuelson [4]. 

Although it was well-known that the distribution of returns or of logarithmic 
returns approximately followed a stable law, there was a barrier to the appli- 
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cation of these concepts in the financial practice. Indeed, stable distributions 
have non-finite variance, and this leads to many mathematical difficulties (for 
a discussion on this point the reader is referred to chapter 3 of Merton's book 
[4]). Therefore, in mainstream finance, both theoreticians and practitioners 
prefer to use the more tractable continuous Wiener process instead of discon- 
tinuous Levy flights. A way of overcoming these difficulties has been provided 
by empirical studies suggesting the use of truncated Levy flights, characterized 
by probability density distributions with finite moments [5-8]. 

In financial markets, not only prices and returns can be considered as random 
variables, but also the waiting time between two transactions varies randomly. 
So far, a large part of the financial practice is based on daily price changes. 
However, a company specialized in intra-day transactions and high-frequency 
data analysis, Olsen & Associates, has published various working papers re- 
lated to the time behaviour of tick-by-tick data (see, for instance, on fractional 
time ref. [9] and on mean first passage time ref. [10]). 

The purpose of this paper is to present a rather general phenomenological 
theory of tick-by-tick dynamics in financial markets. Many well-known aspects, 
such as the Levy scaling form of ref. [6], follow as particular cases of the theory. 
The theory fully takes into account the non-Markovian and non-local character 
of financial time series. Predictions on the long-time behaviour of the waiting- 
time probability density arc presented. Finally, a more general scaling form is 
given, based on the solution of the fractional diffusion equation. 

The paper is divided as follows. In Sec. 2, we discuss the relevance of continuous- 
time random walks in finance by explicitly performing a mapping from finan- 
cial data to random walks. In Sec. 3, we present the master equation and we 
show that it reduces to the fractional diffusion equation in the hydrodynamic 
limit (corresponding to a long jump-observation scale and long observation 
times) if some simple scaling assumptions on the jump and waiting-time prob- 
ability densities hold true. Sec. 4 is devoted to the solutions of the fractional 
diffusion equation and their natural scaling properties. Finally, in Sec. 5, we 
point out the main conclusions and outline the direction for future work. 

As a final remark, let us stress that the theory of continuous-time random 
walks is well developed [11,12]. and its relation to the fractional diffusion 
equation and fractional calculus [13] has been recently discussed by various 
authors [14-18]. However, as far as we know, these concepts have not yet been 
applied to finance in the form we present here. 
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2 Continuous-time random walk in finance 



The price dynamics in financial markets can be mapped onto a random walk 
whose properties are studied in continuous, rather than discrete, time [4]. Here, 
we shall perform this mapping, pioneered by Bachelier [1] and fully exploited 
by Samuelson and his school [4], in a rather general way 

As a matter of fact, there are various ways in which to embed a random walk in 
continuous time. Here, we shall base our approach on the so-called continuous- 
time random walk (henceforth abbreviated as CTRW) in which time intervals 
between successive steps are random variables, as discussed by MontroU and 
Weiss [11] 

Let S{t) denote the price of an asset or the value of an index at time t. In a 
real market, prices are fixed when demand and offer meet and a transaction 
occurs. In this case, we say that a trade takes place. In finance, returns rather 
than prices are considered. For this reason, in the following we shall take into 
account the variable x{t) = log S{t), that is the logarithm of the price. Indeed, 
for a small price variation AS* = S'(ti+i) — S{ti), the return r = AS/S{ti) and 
the logarithmic return riog = log[S{ti+i)/ S{ti)] virtually coincide. 

As we mentioned before, in financial markets, not only prices can be modelled 

as random variables, but also waiting times between two consecutive transac- 
tions vary in a stochastic fashion. Therefore, the time series {x{ti)} is char- 
acterised by (p{i,T), the joint probability density of jumps = x{ti+i) — x{ti) 
and of waiting times Tj = ij+i —ti. The joint density satisfies the normalization 
condition / / d^dT(f{^,T) — 1. 

MontroU and Weiss [11] have shown that the Fourier-Laplace transform of 
p{x, t), the probability density function, pdf, of finding the value x of the price 
logarithm (which is the diffusing quantity in our case) at time t, is given by: 

PM= '-^'-'\ ' (1) 

S 1 — (p[K, S) 

where 

p{k,s)= dt dxe-^* + "^^p{x,t), (2) 

Jo J-oo 

and '0(r) = / d^(p{^,r) is the waiting time pdf. 

Let us now consider the situation in which the waiting time and the size of 
the step are independent. In this case the joint density function, can be 
factorized, namely written as the product of a "spatial" part and a temporal 
part: (^(C,t) = \{^)iP{t). Here A({) is the probability for a displacement ^ 
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in each single step (transition probability density). Now, the normalization 
condition for the transition pdf: / d^\{^) = 1 must be added to that for the 
probability density of the waiting time / dTip{T) = 1. 

As a consequence we get: 

^ ' ' s 1- A(«;)V(s) 1- A(«;)V'(s) ' 

where \{k) , the Fourier transform of the transition probability density, is 
usually called the structure function of the random walk and ^(s) = (1 — 
is the Laplace transform of 

/■oo rt 
^{t) = / ^{t') dt'^1- / 7jj{t') dt' . (4) 

Jt Jo 

is the survival probability at the initial point position {to = 0) [15]. 
jQip{t')dt' represents the probability that at least one step is taken in the 
interval (0,t), hence \E'(t) is the probability that the diffusing quantity does 
not change during the time interval of duration t after a jump [18]. 

According to Weiss [19], "^(t) can be viewed as the probability that the du- 
ration of a given interval between successive steps is strictly greater than t 
and is the peculiar function needed to specify the probability of the displace- 
ment at time t* -|- t in a CTRW, where t* is the instant of the last jump. The 
waiting-time pdf is related to ^{t) by the formula: V'(^) — —d^{t)/dt. 

Let us finally remark that, in general, the CTRW is a non-Markovian model 
[19], as at any time one has to know the value of the diffusing quantity as well 
as the time at which the last step took place in order to predict the further 
course of the walk. The non-Markovian property arises because the time of the 
previous step does vary and could be even t = 0, so that the complete history 
of the process must be taken into account at all times. The only Markovian 
version of the CTRW is the one in which the waiting time pdf, iP{t), is a 
negative exponential: 

V'W = ^exp(-r/r), 

where T is the average time between successive steps. Only for this form of 
the density, the probability that a step of the random walk will take place 
in [t , t + dt) is dt/T , as dt 0, independent of the time at which the 
immediately preceding step occurred. This is not true of any other form of 
V^(t). 
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3 Master equations and fractional diffusion 



The master equation governing the probabihty density profile in a CTRW can 
be derived by inverting the Fourier- Laplace transform in eq. (3). Rewriting 
(3) as 

we obtain 



rt r+oo 

p{x,t)^Sa;o'^{t)+ dt'ip{t-t') / dx' X{x-x')p{x',t') 

Jo J-oo 



(5) 



This form of the master equation is quoted, e.g., in Klafter et al. [20] and Hilfer 
and Anton [15]. However, equivalent forms can be found in the literature. The 
following form shows the non-local and non-Markovian character of the CTRW 
[12,18]: 



d_ 



I rt n+oo 

-p{x,t)^ / dt'(t){t-t') -p{x,t')+ / dx'X{x-x')p{x',t') 

t Jo L J-oo 

here, the kernel is defined through its Laplace transform 



(6) 



1-^(5) 



The above equations allow to compute p{x, t) from the knowledge of the jump 
pdf A(^) and of the waiting-time pdf iP{t). In principle, both these quantities 
are empirically accessible from high-frequency market data, even if, recently, 
within the physics community, emphasis has been given to the jump pdf [6]. 

The time-evolution equation for p{x, t) has a remarkable limit, if some scaling 
conditions on the structure function and on the waiting time pdf are verified. 

Let us assume the following scaling behaviour in the hydrodynamic limit (long- 
jump scale and long observation times): 

A(/«) - 1 - [«[", K^Q, Q<a<2, (7) 

and 

^p{s) s^O, 0<(3<1. (8) 



The above approximations are consistent with the following explicit expres- 
sions for the Fourier and Laplace transforms: 

A(k) = exp(-[K['*) , 0<q;<2, (9) 
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and 



We note that eq. (9) represents the characteristic function for the symmetric 
Levy stable pdf of index a ; for < a < 2 the pdf decays hke as 
|x| ^ oo , for a = 2 the Gaussian pdf is recovered. 

From eq. (10), we observe that 

so that the survival probability turns out to be 

^{t) ^ Ep{-t^) , 0</3<l, (12) 

where 

oo j./3n 

is the Mittag-Leffler function of order /3 [21,22]. Thus the pdf for the waiting 
time is 

^^^^ = "1^^^^ = "1^^^"^'^ ' 0<^^1' (13) 

which is in agreement with the expression obtained in [15] in terms of the 
generalized Mittag-Leffler function in two parameters. 

For < /9 < 1 the Mittag-Leffler function Eg{~f/) is known to be, for t > 0, 
a completely monotonic function of t, decreasing from 1 (at t = 0) to like 
as i — > oo [22] . As a consequence the pdf for the waiting time is strictly 
positive and monotonically decreasing to zero like r(^+^) . For /3 = 1 the 
Mittag-Leffler function reduces to exp(— i) and we recover from eqs. (11 - 13) 
the Markovian CTRW. 

If we insert eqs. (7) and (8) into eq. (1), we get the hmiting relation: 

S^^{k, s) + \K\y{K, S) = . (14) 



Inverting eq. (14), we obtain the timc-cvolution equation for p{x,t) in the 
hydrodynamic limit. If < /3 < 1 and < o; < 2, we have, for x & R: 

dtP ~ d\x\- +r(i-/3)^(^)' (^^")- ^^^> 
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In cq. (15), we have introduced the fractional derivatives /dt^ and /d\x\'^ 
defined as the inverse Laplace and Fourier transforms of and — re- 
spectively [13,17]. Fractional derivatives are non-local operators belonging to 
the larger class of pseudo-differential operators [23,24], which allow power- law 
effects. In particular, the "time" operator in eq. (15) is the Riemann-Liouville 
fractional derivative of order /3 defined as (if < /3 < 1): 

-Jit) 




whereas the "jump" operator is the Ricsz fractional derivative of order a 
which, if < a < 2 can be represented as [13]: 

''°-/M = r(i + c)"°'°-^^' r/(^ + 0-2/W + /(^-0^, 



d\x\" TT Jo ^ 



Finally, let us mention that eq. (14) was derived by Weiss [19] and by Afanas'ev 
and co-workers [25]. Moreover, the above derivation of eq. (15) was implicit in 
a paper by Fogedby [14] and was explicitly presented by Compte [16] and by 
Saichev and Zaslavsky [17]. 



4 Levy flights and scaling of solutions 



We start this section with the analysis of a particular case of eq. (15), the 
limit /? ^ 1, where we have (in the weak sense) [17]: 

lim — — — S(t), 

/3-ir(i-/3) ^ ^' 

and eq. (15) becomes equivalent to the following initial value problem: 

dp{x, t) d°'p{x, t) 



dt ' 



p{x,0) = S{x). (16) 



The Cauchy problem (16) can be solved by Fourier-transforming both sides 
of the equation with respect to x. After integrating and inverse Fourier-trans- 
forming, one gets: 

pM = 1^lJ-^], (17) 



where La{u) is the Levy standardized probability density function: 



1 f + OO 

L^{u) = — / e-'^^'-^^^-dq. (18) 
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Taking the limit /? ^ 1 in cq. (15) corresponds to considering independent 
time increments. Continuous-time random walks whose pdf p{x, t) is given by 
eq. (17) are called symmetric Levy flights or better symmetric a-stable Levy 
processes [3,26]. In 1963, analysing the scaling properties of financial time- 
series, Mandelbrot [2] found that the empirical pdf p{x, t) could be well fitted 
by the Levy density function (17) with o; = 1.7. As we mentioned in the in- 
troduction, the main difficulty in dealing with the Levy distribution is that 
its moments diverge. For < o; < 2, the only bounded finite moments have 
index 7 satisfying — 1 < 7 < a. For this reason, the results of Mandelbrot 
were well-known but not much used in mainstream quantitative finance [4]. 
The recent empirical analysis of Mantegna and Stanley [6] suggests that trun- 
cated Levy flights should be used instead, as good models for flnancial price 
dynamics [5]. Koponen [7] introduced a class of truncated Levy flights, which 
was successively generalized by Boyarchenko and Levendorskii [8]. However, 
all these studies somehow neglected the waiting-time pdf. 

In the general case, the Cauchy problem of eq. (15) can be solved by the same 
technique used above. There is, however, a mathematical subtlety. In order to 
give a meaning to the Cauchy problem, the Riemann-Liouville operator must 
be replaced by the Caputo fractional derivative of order /3 [22,27]: 

^ 1 d ( ft fir) , 1 t-^ 

Now, the solution is: 

- J^W^,P {^) ^ (19) 

and Wa,i3{u) is the following scaling function: 

1 f+00 

WaA^) = TT / e-^"£;^(-|gr)c^g, (20) 

where Ep the Mittag-Leffler function of order /3 and argument z — — \q\'^ 

Further empirical studies on high-frequency flnancial data may reveal the scal- 
ing form (20), if the waiting-time pdf satisfles the asymptotics (8). 



5 Conclusions and outlook 

In this paper, we argued that the continuous-time random walk (CTRW) is a 
good phenomenological model for high-frequency price dynamics in flnancial 
markets, as, in general, this dynamics is non-Markovian and/or non-local. 
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CTRW naturally leads to the so-called fractional diffusion equation in the 
hydrodynamic limit if some scaling properties of the waiting time pdf iP{t) 
and of the jump pdf A(^) hold true in that limit. This point needs a further 
discussion. Indeed, the scaling regime of eqs. (8) and (7) breaks down for very 
large jumps. For this reason, truncated Levy flights have been introduced [5- 
8] . Preliminary investigations [28] on high frequency financial data show that 
a similar problem is present for the waiting time pdf. Nevertheless, we can 
view the fractional diffusion equation (15) as a model for approximating the 
true behaviour of returns in financial markets. 

In the region where the dynamics is well approximated by eq. (15), we expect 
the following scaling for the waiting-time pdf (see the discussion in Sec. 3): 

V'(t) ~ T-^ (21) 

where = P + 1 varies in the range 1 < < 2. Consequently, the more 
complex scahng form (20) should hold true. 

Empirical analyses on market high-frequency data will be necessary in order 
to verify these predictions. In any case, we expect that the concepts of CTRW 
and of fractional calculus will be of help in practical applications such as 
option pricing, as they provide an intuitive background for dealing with non- 
Markovian and non-local random processes. 

In this paper, the mathematical apparatus has been kept to a minimum, the 
interested reader will find full mathematical details in a forthcoming paper 
[29]. 
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